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ABSTRACT 

The  spanwise  variation  in  pressure  due  to  gravity  (occurring 
in  experimental  work)  was  included  in  the  numerical  calculations 
of  a supercavitating  hydrofoil  of  finite  spein.  The  effect  was 
negligible  on  the  lift  and  moment  coefficients  but  a small 
change  in  the  cavity  length  was  noticed  at  the  foil  centroid. 

The  major  change  was  in  the  cavity  shape.  The  cavity  predicted 
by  the  new  numerical  model  was  longer  near  the  root  and  shorter 
near  the  tip  as  was  expected. 
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I.  INTRODUCTION 

To  insure  the  accuracy  of  numerical  calculations,  they 

must  be  compared  to  experimental  results.  In  the  case  of  the 

three-dimensional  supercavitating  hydrofoil  a numerical  model 

has  been  developed  and  has  compared  favorably  to  experimental 

results  (Jiang  and  Leehey  [1]) . However,  the  experimental 

results  were  obtained  using  a technique  in  which  the  foil  was 

mounted  vertically  in  a water  tunnel  with  the  tip  pointed 

down  (see  Fig.  1) . This  causes  a spanwise  variation  of  pressure 

and  therefore  a spanwise  variation  of  cavitation  nximber 
1 2 

(a  = “PjJ/^pUoo  ).  The  experimental  procedure  most  widely 
used  was  to  calculate  the  cavitation  number  based  on  the  free 
stream  pressure  at  the  foil  centroid  and  the  measured  cavity 
pressure  (sometimes  the  vapor  pressure  is  used  as  the  cavity 
pressure) . It  is  the  intent  of  this  paper  to  test  the  effect 
of  altering  the  numerical  calculation  to  include  these  effects. 


II.  THEORY 

A numerical  lifting  surface  model  for  supercavitating 
hydrofoils  of  finite  span  has  been  developed.  It  is  a linear 
theory  for  use  on  thin  wings  at  small  angles  of  attack.  The 
model  uses  discrete  vortices  and  sources  located  on  a lattice 
to  represent  the  foil.  Each  element  of  the  lattice  contains 
a discrete  bound  vortex,  a trailing  vortex  and  a source.  The 
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bound  vortex  is  located  at  the  quarter  chord  line  of  the  element 
and  the  source  is  a constant  distribution  along  the  three- 
quarter  chord  except  at  the  leading  edge  where  the  source  is 
located  at  the  quarter  chord  line.  The  lattice  used  in  this 
model  is  one  obtained  by  dividing  the  foil  and  cavity  spanwise 
into  strips  with  cosine  spacing  and  chordwise  into  strips  of 
constant  spacing  (see  Fig.  2) . The  solution  of  the  problem  is 
obtained  by  reducing  the  coupled  integral  equations  to  a set 
of  simultaneous  algebraic  equations.  The  cavity  lengths  were 
iterated  at  each  chordwise  strip  to  obtain  the  proper  cavitation 
number  over  the  cavity  planform. 

The  foil  and  cavity  surfaces  are  collapsed  onto  the  X-Z 
plane  (see  Fig.  3) . The  jumps  in  tangential  and  normal  per- 
tabation  velocity  components  across  the  foil  and  cavity  are: 

1.  u(x,z,+0)  - u(x,z,-0)  = - y(x,z) 

2.  v(x,z,+0)  - v(x,z,-0)  = q(x,z) 

where  y is  the  vortex  distribution  and  q is  the  source 
distribution.  The  boundary  conditions  are: 

3.  v(x,z,-0)  = dy (x, z ,-0) /dx 

4.  u = cr/2 

5.  q(^,z)d^  = 0 
l(z) 
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where  the  integration  is  done  along  each  strip  of  elements 
from  the  leading  edge  to  the  end  of  the  cavity.  The  induced 
velocities  are  calculated  at  the  midspan,  three-quarter 
chord  position  of  each  element  in  the  model.  The  condition 
of  fixed  cavity  pressure  is  also  tested  in  each  element,  but 
at  the  midspan,  quarter  chord  location.  The  integral  equations 
are : 


\ 


with  the  boundary  conditions  they  become  the  following  set 
of  algebraic  equations. 
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9. 

l/47r  Z 

i/ j 

^jkt  = -1 

i=l,2, 

j=l,2. 

. . .Ml 
. . .N 

10. 

1/217  E b.  ..  .q.  ./a  ~ o/a  = 0 

ij 

i=l,2, 

j=l,2. 

. . .M 
. . .N 

for  k>Ml 

11. 

E q. ./a 
i 

AS.  . = 0 

ID 

j=l,2. 

. . .N 

where  and  sre  the  v and  u velocity  components  at 

the  control  point  (k,£.),  caused  by  a unit  strength  Vortex 
and  source  at  element  (i,j).  M is  the  number  of  elements 
along  the  chord  in  the  cavity  and  Ml  and  N are  the  number  of 
elements  on  the  foil  in  the  chordwise  and  spanwise  directions, 
respectively.  The  accuracy  of  this  model  has  been  tested  by 
Jiang  and  Leehey  [l]  and  results  compare  well  with  both 
analytical  and  experimental  results  (see  also  figures  4-6). 

To  include  the  spanwise  effects  of  gravity  for  the  foil 
in  the  verticle  (experimental)  position,  equations  (4)  and  (10) 
above  must  be  changed  to: 

4 ' . u = o (z) /2 

10'.  + l/2n  1 = 0 i=l,2,...M 

i,  j j=l,2, . . .N 

Yj^£=0  for  lc>Ml 

where  is  now  a function  of  span  position  (S,)  . The  change 
in  free  stream  pressure  AP^  due  to  a change  in  vertical 
position  h is  given  by: 
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12.  AP^  = pgh  (see  Fig.  1). 

Therefore  the  change  in  cavitation  number,  Aa  due  to  a change 
in  position,  h is  given  by: 

13.  Aa  = (P^--P^)/|  pu^2_(p^_p^)/l 

= AP„/|pU„2  = 2gh/U„^ 

Experimental  procedures  quote  cavitation  numbers  based  on  free 
stream  pressures  measured  at  the  centroid  of  the  foil.  There- 
fore the  variation  of  cavitation  number  over  the  span  can  be 
determined  using  equation  (13)  and  the  "overall"  cavitation 
number  (at  the  centroid  of  the  elliptical  planform) . 

III.  RESULTS 

Equation  (13)  was  added  to  the  numerical  technique  and 
the  input  cavitation  number  was  taken  as  the  corresponding 
to  the  foil  centroid.  Then,  an  array  of  cavitation  numbers, 
each  corresponding  to  the  center  of  a strip  of  elements  along 
the  chord  of  the  foil  were  determined  (see  Fig.  2) . The  model 
was  also  changed  so  that  equation  (10')  was  used  in  place  of 
equation  (10) . Several  test  conditions  were  calculated  and 
the  results  are  plotted  in  Figures  4-7.  The  test  conditions 
used  were  for  two  foils  of  aspect  ratio  five  and  elliptical 
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planform,  each  at  four  different  values  of  cavitation  number 
divided  by  angle  of  attack  (a/a) . One  foil  had  a semi-span 
of  10  inches,  the  other  was  15  inches.  (Experimental  results 
are  available  for  foils  of  these  sizes  from  Maixner  [3],  these 
are  shown  for  comparison.  Analytical  results  of  Leehey  [4J 
are  also  shown.)  The  values  of  a/a  used  were  1.0,  1.5,  2.0 

and  2.5.  The  results  show  that  the  spanwise  effect  of  gravity 

is  negligible  in  this  model.  Both  lift  coefficient  over 
angle  of  attack  (C-^/a)  (Fig.  4)  and  moment  coefficient  over 

angle  of  attack  (Cj^^/a)  (Fig.  5)  increased  in  magnitude  over 

the  range  of  a/a  tested,  but  the  increase  was  only  in  the 
third  significant  figure.  The  larger  foil  had  larger  values 
of  Cj^/a  and  Cj^^/a,  but  again  the  differences  were  small.  The 
cavity  length  to  chord  ratio  {i/c)  (see  Fig.  6)  was  larger 
than  the  original  numerical  model  by  about  0.2  for  high  a/a 
and  shorter  by  about  0.2  for  low  a/a,  but  both  models 
agree  well  with  experimental  data.  Figure  7 shows  the 
major  difference  between  the  two  models,  the  cavity  shape. 

The  cavity  shape  is  plotted  for  the  small  foil  at  a/a  = 1.5 
as  it  was  calculated  by  the  two  numerical  methods.  The 
new  model  obtains  a shorter  cavity  near  the  tip  due  to  the 
higher  cavitation  number  there  and  a longer  cavity  near  the 
root  due  to  the  lower  cavitation  number.  The  experimental 
result  for  this  condition  (from  Maixner  [3])  is  also  shown. 
Figure  8 is  the  photograph  used  to  determine  this  result. 
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The  experimental  result  compares  favorably  with  the  new  model 
near  the  tip  but  the  cavity  is  shorter  near  the  root  due  to 
the  boundary  layer  on  the  upper  wall  of  the  test  section. 


IV.  CONCLUSIONS 

It  is  clear  from  the  results  presented  here  that  the 
spanwise  effect  of  gravity  can  be  neglected  for  the  numerical 
model.  The  effect  of  including  the  spcinwise  variation  in 
cavitation  nvimber  on  lift  and  moment  calculations  has  been 
shown  to  be  only  on  the  order  of  one  percent.  The  cavity 
length  calculations  changed  by  about  five  percent  (that  is 
when  the  cavity  length  is  measured  at  the  span  position  of  the 
centroid) . The  most  striking  difference  found  in  the 
numerical  results  when  the  spanwise  variation  in  cavitation 
number  was  taken  into  account  was  in  the  cavity  shape  predicted. 
This  shows  that  the  reason  for  the  insignificant  differences 
in  the  models  is  that  the  cavities  predicted  are  about  the  Scune 
length  near  the  centroid,  where  the  cavitation  number  was 
measured  in  experimental  work. 
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FIGURE  2 NUMERICAL  MODEL 


FIGURE  3 FOIL  PLANFORM  AND  TYPICAL  VORTEX  SOURCE 
ELEMENT 


Moment  Coefficient  of  Elliptical  Planform,  Aspect  Ratio 


FIGURE  6 Cavity  Length  to  Chord  Length  Ratio  of  Elliptical  Planf orm 
Aspect  Ratio  5. 
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FIGURE  7 CAVITY  SHAPES,  ASPECT  RATIO  5.  cr/a 
SEMI -SPAN  * 10"  SCALE  = 1/2  SIZE 


(a)  SMALL  foil:  a*  12*,  Gr*  .3545,  L/c*l.3 


FIGURE  8 Photograph  of  Cavity  Shape. 
(From  Maixner  [3]) . 
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